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We study theoretically the magnetoresistance pxx{B) of a two-dimensional electron gas scattered 
by a random ensemble of impenetrable discs in the presence of a long-range correlated random 
potential. We believe that this model describes a high-mobility semiconductor heterostructure with 
a random array of antidots. We show that the interplay of scattering by the two types of disorder 
generates new behavior of Pxx{B) which is absent for only one kind of disorder. We demonstrate 
that even a weak long-range disorder becomes important with increasing B. In particular, although 
pxx{B) vanishes in the limit of large B when only one type of disorder is present, we show that 
it keeps growing with increasing B in the antidot array in the presence of smooth disorder. The 
reversal of the behavior of pxx{B) is due to a mutual destruction of the quasiclassical localization 
induced by a strong magnetic field: specifically, the adiabatic localization in the long-range Gaussian 
disorder is washed out by the scattering on hard discs, whereas the adiabatic drift and related 
percolation of cyclotron orbits destroys the localization in the dilute system of hard discs. For 
intermediate magnetic fields in a dilute antidot array, we show the existence of a strong negative 
magnetoresistance, which leads to a nonmonotonic dependence of pxx{B). 



I. INTRODUCTION 



In recent years, there has been a revival of interest in 
quasiclassical transport properties of a two-dimensional 
electron gas (2DEG). This has been largely motivated 
by the experimental progress in controlled preparation 
of nanostructured semiconductor systems |Q] and, in par- 
ticular, by the experimental and practical importance of 
high-mobility heterostructures, in which impurities are 
separated from the 2DEG by a wide spacer. On the the- 
oretical side, much of the recent interest in quasiclassics 
on the nanometer scale has been related to the realiza- 
tion that the classical dynamics in a disordered system is 
in fact much richer than the idealized Drude picture sug- 
gests. Indeed, as far as ballistic mesoscopic systems are 
concerned, electron transport has been studied in terms 
of quasiclassical dynamics in great detail However, 
in diffusive systems with smooth disorder, a quasiclas- 
sical treatment of electron kinetics is also appropriate 
and has been shown to lead to new transport regimes. 
To describe the transport properties of such system, one 
sometimes has to completely abandon theories based on 
the Boltzmann equation. In Boltzmann transport theory, 
formulated in terms of a set of relaxation times, quasiclas- 
sics leads to the Drude results: analytical behavior of the 
ac conductivity cr(uj) at w — > 0, zero magnetoresistance 
(MR), etc. It has been demonstrated, however, that qua- 
siclassical memory effects, neglected in the conventional 
Boltzmann approach, yield a wealth of anomalous trans- 
port properties of a 2DEG subject to long-range disorder. 
In particular, non-Markovian kinetics gives rise to a qua- 
siclassical zero-frequency anomaly (see and references 
therein) in the ac response of a disordered 2DEG, asso- 
ciated with return processes in the presence of smooth 



inhomogeneities. Specifically, the return-induced correc- 
tion to Kea{Lo) exhibits a kink cx \lu\. Another manifes- 
tation of non-Markovian kinetics is a strong positive MR 
in low magnetic fields, which is able to explain jH the 
otherwise puzzling positive MR observed near half-filling 
of the lowest Landau level in the fractional quantum Hall 
regime. The strength of the above anomalies depends on 
the ratio d/l, where d is the correlation radius of disor- 
der, I the mean free path, and grows with increasing d/l 
as a power of this parameter. Since quantum corrections 
are governed by a different small parameter l/kpl <^ 1, 
where kp is the Fermi wavevector, it is the long-range 
correlations of disorder with kpd ^ 1 that reveal the 
quasiclassical anomalies. The condition kpd 3> 1 is typ- 
ically well satisfied in high-mobility semiconductor het- 
erostructures. 

In this paper, we consider the quasiclassical magne- 
totransport properties of a 2DEG in a random array of 
antidots (AD). The transport (dc and far- infrared) prop- 
erties of AD arrays, both periodic and random, have been 
the subject of many recent experiments, see, e.g., 
and references therein. In periodic arrays (for a review 
interest has been focused on geometric reso- 
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nances which are associated with the periodicity and re- 
sult, in particular, in commensurability peaks in the MR 
On the other hand, random arrays (see, e.g., 
constitute a remarkable disordered sys- 
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tem where the ADs play the role of hard-wall scatterers. 
We aim to study the MR in random AD arrays and there- 
fore assume that there exist two types of disorder: ADs, 
which we model as impenetrable hard discs that scatter 
electrons, and a smooth random potential, created in the 
heterostructures by charged impurities behind a spacer. 
The quasiclassical MR in each of the limits, where only 
one type of disorder is present, is well understood by now 
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(see Sec. g). The purpose of the paper is to demonstrate 
that the interplay of the two types of disorder yields new 
physics that is absent in the limiting cases. We will show 
that, although in the extreme of a strong magnetic field 
B —^ oo the dissipative resistivity Pxx{B) tends to zero 
in either of the limiting cases, it diverges in the presence 
of both types of disorder. In particular, in the experi- 
mentally relevant situation of relatively weak long-range 
disorder, i.e., when the mean free path at zero B is de- 
termined by scattering on ADs, the presence of the weak 
long-range fluctuations will nonetheless become of crucial 
importance with increasing B. It is worth noting that our 
model can also be applicable to the description of the MR 
in an unstructured 2DEG with residual interface impu- 
rities playing the role of antidots (large-angle scattering 
on residual impurities is known to become important in 
unstructured samples with a wide spacer 

The paper is organized as follows. We give a brief re- 
view of past work on the quasiclassical MR in Sec, p^. In 



the body of the paper, we first consider in Sec. Ill the 



MR at a moderately strong B, when the collision time 
for scattering on ADs is not affected by the magnetic 
field. Then, in Sec. IV, we turn to the limit of strong B, 



where the collision time is renormalized as compared to 
the Drude value and, in the extreme of very large B, a 
single act of scattering involves "skipping" of cyclotron 
orbits along the surface of ADs. The whole picture turns 
out to be rather complex and we choose the following 
logic of presentation. We fix the zero-B mean free paths 
for scattering on ADs and on the long-range disorder and 
for different values of the density n of ADs sweep the 
magnetic field. In Sec. [II we start with the "hydrody- 
namic limit" of infinite n and then gradually decrease 
n. In Sec. IV we first consider a single act of scatter- 



ing on an AD for large B, then proceed to analyze the 
strong-B transport in an AD array. We present results of 
numerical simulations in Sec. |^ and summarize in Sec. VI 
[where the qualitative behavior of Pxx{B) is illustrated in 
Fig.H]. 



II. 



OUTLINE OF KNOWN RESULTS: LIMITING 
CASES 



A. Lorentz model 



Is/Rc, where Rc is the Larmor radius, and reads ||2| 
units of the zero-B resistivity po, 



PxxjB) fh_ 
Po \Rc 



(1) 



with F{0) = 1 and Fix > 1) ~ 97r/8a; < 1. In Drude 
theory, the dissipative resistivity is not affected by a mag- 
netic field and F{x) — 1 for all x. The nontrivial kinetic 
problem (|l|) is in fact fully solvable and the exact expres- 
sion for the conductivity tensor at arbitrary 1$ / Rc can be 
found in |25|j2^ ] (see also ||2^,^ for numerical simulations 
of the problem). 

The falloff of pxx B^^ is related to the peculiarity 
of the Lorentz model: at finite B, there are electrons 
that move freely in steady cyclotron orbits and never 
hit a scatterer; those electrons do not contribute to pxx 
and their density grows with increasing B. The conduc- 
tivity is then due to electrons that experience multiple 
collisions with a scatterer by moving in "rosette" orbits 
around it (Fig. |l|) until they hit another scatterer, which 
results in a diffusive hopping of the "rosette states" . At 
finite concentration n, the Lorentz model has a metal- 
insulator transition |^ , ^ at Rc ^ n^^l"^: for larger B 
the dissipative conductivity is strictly zero, as shown in 
Fig. 0. 




FIG. 1. Schematic behavior of the magnetoresistivity 
Pxx{B) as a function of Is/ttRc in the Lorentz model. In- 
set: Rosette orbit of an electron bound to a hard disc (shown 
by the shaded circle) in a magnetic field. 



We start by briefiy recalling the known results for the 
classical Lorentz model in two dimensions (hard discs of 
radius a, randomly placed with a concentration n; we 
assume that na^ ^ 1 and kpa ^ 1, so that the mean 
free path Is = 3/8na). This is a good model for an AD 
array in a heterostructure. As was pointed out in p^ , p5t , 
Drude theory fails completely to describe magnetotrans- 
port in this system. In the limit n — s- oo, a — > 0, with 
Is held fixed, the resistivity depends on a single variable 



B. Long-range disorder 

Now let us recall what is known about MR in the 
case of a smooth (allowing for a quasiclassical treatment) 
Gaussian (in the sense of statistics of fluctuations) ran- 
dom scalar potential. There are two sources of quasiclas- 
sical MR (we consider elastic scattering on an isotropic 
Fermi surface). 
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First, note that the MR is strictly zero in Boltz- 
mann theory only in the limit of white-noise disorder, 
whereas if disorder is correlated on a finite spatial scale 
d, the collision-integral approximation allows for a fi- 
nite MR 1^,^, due to a cyclotron bending of trajec- 
tories within this correlation radius. This simple effect 
is governed by the parameter d/Rc [at small B it yields 




FIG. 2. Schematic behavior of the magnetoresistivity 
pxx{B) as a function of d/Rc for a Gaussian smooth random 
potential. 

Second, there is MR associated with memory effects 
and, to calculate this, one has to go beyond the collision- 
integral approximation. The memory effects are brought 
about by correlations of scattering acts at the points 
where quasiclassical trajectories self-intersect. These ef- 
fects give the main contribution to llp^x/ Po at large 
enough B, where the governing parameter \s d/5 with 
5 being a characteristic shift of the center of a cyclotron 
orbit after one revolution. For Rc/d > 1 the shift is (see 
pO| and references therein) 



Rc{Rc/Il 



,1/2 



Rc>d 



(2) 



where II is the mean free path in the smooth random po- 
tential (experimentally, iLjd ^ 10^ — 10'^ in high- mobility 
samples). According to lS.pxx/ Po ^ < 1. The 

return-induced contribution becomes much larger than 
that related to the effect of B on the collision integral 
at {5/Rcf < d/5, i.e., at Rc < iLid/hf^^. The exact 
expression for the MR in the hmit {S/Rc)^ ^ d/(5 <C 1 
in a hetcrostructure with a spacer d is 



Ap,4B) ^ 2C(3/2) 



9/2 



(3) 



This equation is valid with increasing magnetic field up 
to d/S ~ 1, where Ap^xiB)/ po becomes of order unity. 
At higher fields, when 6/d <^ 1, the strong positive MR 



is followed by a sharp (exponential) falloff of pxx with 
growing B [R0| : 



In 



Pxx 

Po 



2/3 



(4) 



which is due to the increasing adiabaticity of the elec- 
tron dynamics and the related quasiclassical localiza- 
tion po| ,p| . The self- intersection induced MR, given by 
Eq. (|^) , may be considered as a precursor of the adiabatic 
localization. In the limit of large B, when Rc/d becomes 
small, d is given by 



Rl/idh 



xl/2 



Rr <d 



(5) 



The nonmonotonic behavior (||),(^ of the MR in the case 
of a purely Gaussian long-range random potential is il- 
lustrated in Fig. ^. 

To conclude the brief overview, it is worth noting that 
in the limit of weak inhomogeneities the return-induced 
MR depends in an essential way on the behavior of the 
disorder under time reversal. In particular, it is strongly 
enhanced in the case of a random magnetic field M . 



III. EFFECT OF CYCLOTRON DRIFT ON 
TRANSPORT IN ANTIDOT ARRAYS 



A. Parameters of the problem 

We now turn to the MR in the presence of both a 
random array of hard discs and long-range Gaussian dis- 
order, which we characterize by the mean free paths at 
zero magnetic field Is and Z^, respectively. We assume 
that Is/Il ^ 1, which describes a typical experimental 
situation. As in Sec. ||, n will denote the concentration 
of ADs, a their radius, d the correlation length of the 
smooth random potential, Rc the cyclotron radius, and 
6 the characteristic shift, due to scattering on the long- 
range disorder, of the cyclotron orbit after one revolution. 
Throughout the paper we assume d/a^ 1. 

We are interested in strong effects in the behavior of 
Pxx{B): for Is/Il ^ 1, these can only occur if Rc/ls ^ 1- 
Moreover, for the most part of the paper (namely, with 
the exception of Sec. IV A ), we consider magnetic fields 
which are sufficiently strong in the sense that 6/d < 1. In 
this case, the motion of electrons is characterized by rapid 
cyclotron rotation around the guiding center and slow 
drift of the latter along equipotential lines of a smooth 
random potential. Most of these lines are closed, which 
leads to localization of particles trapped on them. The 
effect of scattering by ADs is to induce transitions be- 
tween equipotential contours and, in this way, allow the 
localized particles to escape. 
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B. Hydrodynamic limit 

Let us first consider the "hydrodynamic hmit" 
(n oo, a 0, Is — const). Clearly, in this limit, 
the effects yielding the falloff of p^x oc B^^ [Eq. (|l])] are 
washed out by infinitesimally weak long-range disorder. 
One might think that then the Drude formula works and 
Pxx{B) / po — 1 for all B. In fact, however, this is not true 
and even a small {Is I^l <C 1) amount of smooth disorder 
becomes a relevant perturbation with increasing B. In- 
deed, in the limit of large B (namely, for 6/d <C 1), the 
problem can be mapped onto that of advection-diffusion 
transport i.e., of a Brownian motion with a diffusion 
coefficient Dq in a spatially random velocity field v(r) 
("steady flow") with V • v = ("incompressible fiuid"). 
In this mapping, the field v(r) describes the adiabatic 
drift of guiding centers of cyclotron orbits due to long- 
range inhomogeneities and Do ~ Rc/'''S} where rg is the 
momentum relaxation time for scattering on ADs. The 
result for the effective (macroscopic) diffusion coefficient 
D in the advection-diffusion problem ||3l| is 



(6) 



if Vdd > Do and D = Dq otherwise. Here Vd is a 
characteristic amplitude of the fiuctuations of v(r) [see 
Eqs. (|Tl|),(p^ below]. Hence the conductivity will be 
strongly enhanced by even a weak long-range disorder 
provided Vdd/Do 3> 1. Since this parameter is a grow- 
ing function of B, the effect of smooth disorder is am- 
plified by the magnetic field. The reason is percolation 
of cyclotron orbits through long-range inhomogeneities: 
the percolation-dominated D can be written as a product 
VdW, where w dis a characteristic width of links of the 
percolation network. The equation w ~ d{Do/vdd)^^^^ in 
the advection-diffusion problem comes from the condition 
of connectivity of the network w^Vd/ L[w) ~ Dq, where 



L{w) - did/wy/^ 



(7) 



is a typical length of the network link [Q. Note that 
the size ^(w) of the elementary cell of the percolation 
network [i.e., a characteristic end-to-end distance for the 
link of length L{'w)] scales as 



^{w) - d{d/w) 



4/3 



(8) 



Although the advection-diffusion model has become 
popular for the description of transport in the high-_B 
limit (in particular, in the quantum Hall regime [ ^ , ^ ), 
we should be careful to check if the scattering on ADs 
can actually be described in this model in terms of the 
diffusion coefficient Dq. Clearly, this requires that w be 
larger than a hopping length for the diffusion process, 
which means w S> -Rc- While this condition is satisfied in 
the extreme of large B, a nontrivial transport regime may 
occur with increasing B in which D 3> Dq but w -C -Rc 



(as we will see below, this is the case if the long-range 
disorder is not too weak) . In this regime, the main contri- 
bution to D comes from electrons that move freely along 
the critical links: the "ballistic" motion along the per- 
colating path is contrasted with the transverse (across 
the drift trajectory) diffusion in the advection-diffusion 
regime. In other words, the number of collisions with 
ADs during the drift along a critical link of the percola- 
tion network is now of order unity. As in the advection- 
diffusion regime, the number of passages of the network 
link between two consecutive changes of critical cells is 
also of order unity. It follows that w obeys the simple 
scaling L{w) ^ VdTs, so that the result for D ^ VdW is 



D ^ Vdd{d/vdTsf'^ 



(9) 



which should be compared, as in the advection-diffusion 
problem, with Dq: Eq. (^ is valid when D > Dq. Note 
that in this new regime D does not contain the hopping 
length, which may be even larger than d. 

We are now prepared to calculate the MR. The chaotic 
scattering on the long-range potential crosses over into 
the adiabatic drift with increasing B at -Rc of order 



i?c = dih/dY'^ , 



(10) 



where 5/d becomes of order unity [cf. Eqs. (||),(^)]. At 
this field, the MR is still weak and transport is completely 
determined by scattering on ADs, whereas at larger B we 
can already use the high- field formulas (^,(||) and write 
Pxx{B)/pq ~ D/Dq. The characteristic drift velocity Vd 
that should be substituted into Eqs. (||),(^ reads (see. 



e.g. 



-Rc ^ d 

Vd = Vp s 



[dlLY'^ \R. 



(11) 



where vp is the Fermi velocity and the function s{x) is 
given by 



x^/^ ^ x < 1 
1 , a; > 1 



(12) 



Which of Eqs. (^),(^ should be used depends on the ratio 
w/Rc, as explained above. Remarkably, the ratio D/Dq 
for both Eqs. (||) and (||) depends only on two variables, 
X = d/Rc and VdTs/d — px^^s{x), where 



Is 



(13) 



As a result, the behavior of the MR as a function of B 
depends, at a — > 0, on the single parameter p, so that we 
can write 



Po \Rc 



(14) 
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Notice that even though we consider the case Is/Il ^ 1, 
the parameter p may be large since the long-range disor- 
der is weak, i.e., d/lL ^ 1- 

At p ^ 1, the MR remains small, i.e., f{x,p) ~ 1, 
for all X <^ p^^l^ . On the upper boundary of this inter- 
val, the percolation starts to renormalize D in accordance 
with Eq. (|^), which yields a power-law growth of PxxiB) 
with further increasing B: 

f{x,p) ^ , p-1/3 « ^ « 1 ; (15) 

^ 1 « a; « p3/io . (16) 

The scaling behavior changes between Eqs. (|l5|) and ( p^ ) 
at a; ~ 1 because of the change in the dependence of Vd 
on B at Rc/d ^ 1. At still larger B, D obeys Eq. (|), 
which gives 



/(a;,p)~(px)i°/i3 , p3/io<<^ 



(17) 



Equations (|l5|)-(|l7[) are illustrated in Fig. 




10/13 



P~' d/Rc 

FIG. 3. Schematic behavior of the magnetoresistivity 
pxx(B) as a function of d/Ra in the hydrodynamic limit 
n — * oo, a Q, Is ~ const for a) p = ls/\/dlL 3> 1 and 
b) p < 1. 

At p ^ 1, the range of x where the enhancement of 
the conductivity is described by Eq. (j^) shrinks away, so 
that f{x,p) ~ 1 for all x <C p~^ and behaves according 
to Eq. ( p^ ) at larger x (see Fig. ^). This establishes 



the meaning of the parameter p: if p ^ 1, the Drude 
regime does not match with increasing B the advection- 
diffusion regime (^ directly, but through the interme- 
diate "one-hop" regime (^, whereas if p is large, this 
intermediate phase is absent. Equation (^7|) tells us that 
Pxx{B) (X B^'^/i^ at _B oo. The divergence takes place 
whatever the ratio Is/Il, even if the long-range disor- 
der is weak and does not play a role at i3 = 0. This 
behavior differs drastically from that given by either of 
Eqs. (|1|),(|). 



C. Finite density of antidots 

So far, in Eqs. (p^-dl^), the scattering on ADs has 
been characterized by Is only, through the single param- 
eter p [Eq. (^3|)], which implies the hydrodynamic limit 
n — > oo, a ~f 0. Now we take into account finite-ri ef- 
fects. We begin as before with the case of large p. New 
relevant dimensionless parameters appear, in particular, 
nRcd. Also, since for a fixed ^5 decreasing n means in- 
creasing a, the parameter 5/a may become relevant, in 
which case the scattering on ADs will be affected by the 
magnetic field in an essential way and the collision time 
will not be given by T5 . We will consider effects governed 
by the parameter 5/a in Sec. iV, Until then, let us as- 
sume that (5/a is sufficiently large, so that this parameter 
plays no role for typical electron trajectories. 

Clearly, li 5/a ^ 1, stable rosette-states 



4p5[ are 

still destroyed by the scattering on the long-range poten- 
tial. Naively, one could think that scattering on ADs is 
then chaotic (no trace of the rosette-state dynamics) and 
Eqs. (p^-(p^) apply. In actual fact, provided nRcd <^ 1, 
multiple collisions with the same AD do occur even for 
(5/a ^ 1, as we will see below. At large p, the former con- 
dition is satisfied with increasing B before 5 /a gets small. 
Multiple returns in a dense AD array become possible 
because of the adiabatic localization, which develops at 
5/d<l. 

Let us start by considering the drift regime under the 
condition nR^d ^ 1, Rc/d ^ 1. Typical trajectories of 
guiding centers are closed loops of size ~ d, which means 
that trajectories of electrons circling along cyclotron or- 
bits are bound to within thin rings of width ~ d and 
radius Re- The area of a strip between the inner and 
outer radii of the rings is ~ Red and, if nRcd ^ 1, in 
most rings there are no ADs. Electrons in these rings 
are adiabatically localized and do not, in the adiabatic 
approximation, contribute to Pxx{B) (we will consider 
the possibility of nonadiabatic decay of these states in 
Sec. 
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There are, however, rare rings with a sin- 
gle AD. For electrons in these rings, a typical time ts 
between collisions with ADs is much shorter than ts ■ In- 
deed, the number of cyclotron revolutions before return- 
ing to the region of size 5 around the AD is typically d/5^ 
while the probability of hitting the AD during one such 
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sweep is a/5. It follows that the number of cyclotron 
revolutions before the electron hits the AD is ^ d/a, i.e., 



Ts ~ Rcd/vpa 



(18) 



which gives ts/ts ^ nRcd <C 1. 

Now, a single collision with an AD does not lead to 
a breakaway from the AD. In fact, the electron experi- 
ences multiple collisions with a single AD and each time 
the center of the ring in which the electron is circling 
hops a distance ^ Re', one can visualize this process as 
a random hopping of the center of the ring on a circle 
of radius Rc around the AD (see Fig. |4|). The electron 
"sticks" to the AD for a time much longer than fg. This 
somewhat intricate dynamics of "hopping rings" reminds 
the evolution of the rosette states [ P4|j25[| : in effect it is 
the adiabatic localization in the long-range potential that 
preserves the character of the rosette-state dynamics. A 
breakup will eventually happen when the electron picks 
up a very rare ring containing two ADs (the existence 
of such rings in the area Rc x Rc implies that nR^ ^ 1, 
which we assume in this derivation) . It is straightforward 
to check that the number of scatterings on a given AD 
before getting to another one is ^ 1/nRcd ^ 1. Multi- 
plying the latter by , we find that the time it takes the 
electron to change ADs (separated by a distance ~ Rc) 
is ~ Ts, which yields the diffusion coefficient of electrons 
participating in this type of transport ~ Do, the same 
as in Drude theory. These electrons, however, represent 
a small fraction of the total number of electrons, namely 
^ nRcd. Hence, the contribution to the macroscopic dif- 
fusion coefficient from drift orbits of a characteristic size 
d is ~ DonRcd <C Do- 




FIG. 4. A cartoon picture of scattering of a cyclotron orbit 
on an antidot ata<C<5<Cd<C-Rc^ {nd)~^ . The position of 
the antidot is shown by a cross at the intersection of the rings 
(radius Rc, characteristic width d) which represent the area 
"covered" by the drifting cyclotron orbit. The arrows denote 
hopping of the guiding center of the orbit between drift trajec- 
tories shown by small loops. Provided nR^ ^ 1, the particle 
will typically break away from the antidot when it picks up a 
ring containing one more antidot and hits the latter. 



Having got the contribution of typical trajectories of 
size d, we should take into account that upon hitting an 
AD the particle may hop onto a drift trajectory of size 
^ larger than d. To put it another way, although most 
trajectories with ^ in the interval d <^ ^ {nRc)~^ 
are adiabatically localized, some of them hit ADs and 
mix with the short-scale trajectories of size d considered 
above. This mixing increases the total fraction of de- 
localized trajectories. To calculate the latter, note that 
the probability density P{C) for a point to belong to a 
drift trajectory of size ^ ^ d (we define ^ as a character- 
istic radius of the area to within which the trajectory is 
bounded) scales as 



pio - d/e 



(19) 



No critical exponents are involved here. One way to get 
a quick proof of this is to realize that, according to per- 
colation theory, for zero altitude on a relief map of a 
random landscape, the number of contours of radius ^ 
in the area ^ x ^ is of order unity. It follows that the 
(integrated over the altitude) fraction of space occupied 
by contours of size ~ ^ is ~ Lw/^"^ ^ d/^, where we used 
Eqs. (0),(|), which yields Eq. (|l|). Thus, the fraction 
of trajectories that are delocalized due to collisions with 
ADs is evaluated by integration 



d£, P(0 min{ni?cC, 1} ^ nRcd h\{l/nRcd) , (20) 



which gives merely an additional logarithmic factor. One 
sees that the characteristic ^ < [nRc)^^ are within the 
limits of applicability of the derivation <^ Rc (trajec- 
tories with larger give rise to a percolative contribu- 
considered in Sec. 



tion to pxx considered in Sec. IIIB). Notice that drift 
trajectories that do not hit ADs may be infinitely ex- 
tended only with zero measure and thus do not con- 
tribute to D. Accordingly, D is evaluated as a diffusion 
coefficient of electrons delocalized by the scattering on 
ADs. Since the relevant ^ <^ Rc, the characteristic hop- 
ping length associated with the change of ADs by these 
electrons is Rc. The characteristic rate of hopping be- 
tween two different ADs is given by n (dS/dt) , where 
(dS/dt) - / d^P{^) [A{^)/fsiO] is the average rate at 
which the area explored by the particle stuck to an AD 
grows in time. Here 



fsiO ~ fs{d)ad 



(21) 



[with Ts{d) defined in Eq. (^_8|)] is the time the particle 
resides on a trajectory of size f before being scattered 
out by the same AD, and A{^) ^ i?cC is the area probed 
during this time. These expressions for fs{(,) and A{^) 
are valid for ^ <^ d{S/a)^^^, whereas at larger ^ the par- 
ticle is scattered out before it comes full circle around 
the closed trajectory and both quantities do not depend 
on ^. We see that A{^)/fs{£,) does not depend on ^ and 
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thus (dS/dt) ~ vpa is determined by ^ ^ d. It follows 
that the hopping rate does not change with the inclusion 
of long trajectories and is given by t^^. Accordingly, 
the diffusion coefficient of delocalized particles is ~ Dq. 
We finally get for the macroscopic diffusion coefficient 
D - DonRcd ln(l/ni?crf), or, for the MR: 



PxxjB) 

Po 



nRcd In 



nRrd 



(22) 



We thus see that, despite 5/ a ^ I, the resistivity is 
strongly suppressed as compared to the Drude result at 
nRcd <^ I. We will return to this regime in Sec. fVE, 
where we will show that the actual condition for Eg. (|22|) 
to be valid is {nRcd)^^/'^ , whereas at smaller 5/a, 

I ^ i5/a ^ (nRcd)"^^^ , a slight modification, namely in 
the logarithmic factor in Eq. (p2|), is necessary. 




nd'' 



d/R 



cm 



-IB 

P d/R, 















b) / 




\ ^/ 





nd^ 



(nd'r 



-1/3 



d/Rc 



FIG. 5. Schematic behavior of pxx{B) as a function 
of d/Rc for intermediate magnetic fields at large p and 
{d/lLf'"^ < nd^ < p-i/^ for a) d/i?cm < {nd^f'^ and b) 
d/Rcm S> {nd^Y^^ . The position of the minimum of pxx{B) 
in Fig. |[i is given by d/R^m ~ p-^'^^[nd^ \n{l / nd^ p'^ ' ' . 
Decreasing n leads to the nonmonotonic dependence of 
pxx{B) (cf. Fig. ^ where n oo). The dashed line in Fig. |^ 
shows the behavior of pxx{B) given by Eq. (p3). 



The falloff of the MR described by Eq. (g2|) is illus- 
trated in Fig. ||a. Shown here is also the percolative 
growth of Pxx{B) into which the falloff crosses over at 
sufhciently large B. As is clear from Fig. ||a, in the 
above derivation we implicitly assumed that there ex- 
ists a range of B within which the percolative contri- 
bution is smaller than that given by Eq. (p2[). Indeed, 
since the percolative MR grows with increasing B (see 
Fig. ^), n should be small enough for the B~^\n.B falloff 
not to be masked by the percolation. Let us formu- 
late the condition for the existence of the minimum in 
the dependence of Pxx{B) for p f. If one neglects 
(as everywhere in this section) the nonadiabatic decay 
of drift trajectories, the condition is nd^ <C p~^^^, as 
can be seen from Fig. By matching Eqs. ( p2|) and 
( |l5| ) we find that the minimum occurs at Rc ^ Rem, 
-4/i9[n(i2ln(l/ndV/^)]^/'^- We will 



P 



by taking the nonadiabatic decay 



where d/Rcm 
show in Sec. Ill E , 
into account, that the condition for the existence of the 
minimum actually reads max{nd^, (d/Z^)^/'^} ^ P~^^^, 
which means that Fig. ||a correctly describes the MR if 
{d/hV/^ « nd^. 

To conclude this section, note that the interval of valid- 
ity of Eq. ( ^ ) shrinks to zero if Rc/d < f and, therefore, 
the above considerations describe the behavior of the MR 
only for nd^ ^ 1. According to Fig. ^, for large p the 
nonmonotonic behavior of Pxx{B) develops for nd^ ^ I. 
In the case of small p the picture is similar but the pa- 
rameter S/a becomes relevant, which will be considered 
in Sec. Fv. 



D. "Metal-insulator transition" 

As noted above, the validity of the derivation of 
Eq. (^ ) requires that nR^ ^ I. This condition appears 
already in the Lorentz model: if the opposite inequality is 
satisfied, the system without long-range inhomogeneities 
would be insulating {pxx = 0) ■ In the presence of long- 
range disorder, the percolation mechanism of transport 
prevents pxx from vanishing even at small nR^. How- 
ever, the parameter nR^ determines the interval of B 
where the nontrivial mechanism of diffusion that leads 
to Eq. ( p^ ) is operative, namely I <C nR^ ^ Rc/d. At 
nRc < 1, this mechanism is switched off in a manner 
inherent in a continuous phase transition by formation 
of disconnected clusters of trajectories, very much simi- 
lar to the metal-insulator transition in the Lorentz model 
p^ , p5[ . Hence, the "short-scale" (as opposed to the per- 
colative) MR [Eq. (p^)] behaves near the transition ac- 
cording to 

'Be - B' 



Pxx{B) 



PO 



(„d2)l/2i^ 1 g 

nd"^ 



Be 



(23) 



where G{x) vanishes as a power law at a; —> on the 
conducting side. The critical point B = Be corresponds 
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— 1/2 

to the critical concentration n — Ur ^ Rc ' . This be- 
havior of pxx{B) is shown in Fig. |^. Comparing Figs. |^ 
and ^D, we see that the critical falloff (^3|) is not masked 
by the percolation provided (nd^)^/^ <C d/Rcm- If this 
condition is satisfied, the minimum in the dependence of 
Pxx{B) occurs at nR^ ~ 1. 

We conjecture that the exponent t in the function G{x) 
can be found by mapping the problem of percolation of 
skipping cyclotron orbits onto that of percolation of the 
electric current through an ensemble of conducting cir- 
cles of radius Rc scattered randomly with the density n 
(note that two rosette orbits of radius 2i?c do not mix 
with each other if the distance between the centers of the 
rosettes exceeds 2Rc). The latter problem belongs to the 
universality class of a two-dimensional percolation with 
a finite threshold, for which many critical exponents are 
known (see, e.g., it is worth noting that the percola- 
tion of drift trajectories considered above does not belong 
to this class). In particular, the fraction of space occu- 
pied by the infinite cluster of connected circles vanishes 
near the threshold as (n — with (3 ~ 0.14, whereas 
the conductivity through the infinite cluster exhibits a 
power-law behavior with the critical exponent t ^ 1.2. 



D ^ Di 



vpda 



(25) 



Note that D in this regime is proportional to a product 
of two diffusion coefRcients, Dq and Dna- If, by con- 
trast, fs ^ (P/Dna, the particle will miss this AD and 
will go on exploring phase space in a chaotic way, which 
gives D = Dq. Since (P/Dna at Rc ^ Rc is of the order 
of the cyclotron frequency vp/ Rc, the crossover between 
these two regimes takes place with increasing B when 
Rc <C Rc, (logarithmically) deep in the drift regime. 
Hence, if nRcd <C 1, the nonadiabatic decay of drift tra- 
jectories stretches the region of a chaotic diffusion to the 
point at which DnaTs/d^ ^ 1, which occurs at dj Rc only 
logarithmically larger than {d/lLY^^. At larger S, the 
resistivity starts to fall off sharply, according to Eq. (p^), 
until DnaTs/d^ becomes of order nRcdln{l / nRcd) , where 
this exponential behavior crosses over into the power-law 
falloff described by Eq. (p^. The characteristic features 
in the behavior of Pxx{B) at p ^ 1 associated with the 
nonadiabatic decay of drift trajectories are illustrated in 
Fig. I for the range of B which corresponds to the falloff 
of pxx{B) in Figs. |a,b. At the point d/Rc ^ {d/lLy/^ 
shown in Fig. ^, the adiabatic localization starts to de- 
velop. 



E. Nonadiabatic decay 



In Sec. IIIC, we inferred the B^^ \nB falloff of Pxx{B) 
[Eq. (|2^)] by assuming that the drift picture is applica- 
ble in the whole range nRcd <C 1. This is legitimate if 
nRcd ^ 1, where Rc is defined in Eq. ([l^). Otherwise the 
Drude value of pxx = Pa holds with increasing B up to 
the field where the adiabatic dynamics starts and there is 
an exponentially fast crossover between pxx = Po and pxx 
given by Eq. (22), which is governed by the nonadiabatic 
scattering. 

Let Dna be the diffusion coefficient across drift trajec- 
tories due to their nonadiabatic mixing. Since the rate 
of nonadiabatic transitions depends exponentially on the 
parameter d/S and is therefore locally a wildly fluctuating 
quantity, we should be more specific here: we define Dna 
through a typical time d^ /Dna that it takes to change 
two typical drift trajectories of size ~ d separated by a 
distance ^ d (for thus defined diffusion coefficient In Dna 
scales as d/5, see, e.g., pO|,p|). If nRcd ^ 1, a particle 
which initially resides on a ring with no AD will reach a 
ring containing one in a time 



{d^/ Dna) /nRcd ■ 



(24) 



Then there are two possibilities. If fg <C d^ / Dna [with fs 
defined by Eq. (jl^)], the particle will hit this AD, so that 
collisions with different ADs will occur at a rate ^ t~^, 
which yields a contribution to the macroscopic diffusion 
coefficient D R^/rna ~ D^Dnafs/d'^, i.e.. 




exp (- const, x B ) 



FIG. 6. Schematic representation of the behavior of pxx{B) 
at nd^ < (d/h)^^^ < p~^^^ < 1 for moderately strong 
B which correspond to the falloff of pxx{B) in Figs. ^,b. 
The thick dashed line shows pxx{B) without taking the nona- 
diabatic decay of drift trajectories into account. Provided 
nd^ <^ [d/lh)^^^, the nonadiabatic transitions stretch the 
range of B where pxx[B) ~ po- Beyond this range, pxx{B) 
falls off exponentially with increasing B until this falloff 
crosses over into the _B~^ ln_B behavior. 

Comparing Figs. |^a,b and ^ one can formulate the con- 
dition for the existence of the dip in the dependence of the 
MR on S for p > 1. If n(P > (cJ/^l)^^^, the dip exists in 
the case of nd'^p^^^ ^ 1, as in Figs. ^,b. If, by contrast, 
nd^ <^ (d/lL)^/^, the nonmonotonic behavior in Pxx{B) 
shows up for dp/lL ^ 1. Figure || illustrates the behav- 
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ior of Pxx{B) in the case when the adiabaticity of motion 
in the long-range potential starts with increasing B well 
before the crossover to the percolative growth of Pxx{B). 
This is the case at id/h)^^^ < niin{d/i?cm, (nd2)i/2}. 
It is clear, however, that the above considerations of the 
nonadiabatic decay are equally valid for the case when 
the opposite inequality is satisfied. The only difference 
is that, if (d/Zi)^/^ > niin{(i/i?cm, (j^d^)^/^}, the expo- 
nential falloff of pxx [B) crosses over into the percolative 
growth directly, without passing through the intermedi- 
ate B~^\nB regime. 



IV. STRONG MAGNETIC FIELD: 
RENORMALIZATION OF THE COLLISION 
TIME 

Now let us turn to the extreme of strong B, where the 
parameter 5/a becomes relevant. As will be shown be- 
low, the effect of small 5/a is twofold. First, the time 
t'q between collisions with two different ADs for parti- 
cles drifting along percolative trajectories (which is T5 
at large enough (5/a, how large — see below) gets longer. 
Second, a typical hopping length Rh across the drift tra- 
jectory due to the scattering on an AD (which is Rc at 
1) gets smaller. Therefore, before calculating the 
MR at small 5 /a we have to consider ho w the e lementary 
scattering acts are modified. Belo w, Sec. [V A deals with 
the scattering time r^, Sees. IV B and [V C with the hop- 
ping length Rh- Section [VP studies the MR for small 
6/a. 



^ diver's/ df I' 



(28) 



A. Effective scattering time 



Let us start by considering the case (5/a ^ 1. To find 
Tg, notice that at 5/a <C 1 the length of the drift tra- 
jectory between collisions L5 ~ VdTg cannot depend on 
a, i.e., only the concentration n matters. In the simplest 
case Rc ^ d, the length L5 clearly obeys the equation 



nLsRc ~ 1 
which is rewritten as 

Tg Tsa/5 



(26) 



(27) 



At Rc ^ d, however, there are different regimes for r'g. 
Namely, Eq. (|2^) is valid for r'g at Rc^ d only as long as 
Ls ^ d, i.e., when the drift trajectory between two colli- 
sions can be approximated as a straight line. If L5 3> d, 
which is the case at nRcd <^ 1, the trajectory exhibits 
fractal dimensionality on the scale of Ls. Specifically, the 
length of the percolating trajectory Ls and the distance 
^5 from the starting point are related to each other by 
~ d{Ls/df/^ [see Eqs. (|),(|)], i.e., 



It follows that the cyclotron orbit passes many times 
through the same spatial regions, which increases the 
time Tg. Let first Ls ^ d but ^s Rc [i.e., d ^ 
Ls <C d{Rc/dy^^]. In this regime, S,s [in contrast to 
Ls in Eq. (26)] is of order {nRc)~^: 



which gives 



n^sRc 



a 

TS T 



,5 {nRcdf/^ 



(29) 



(30) 



Now let ^5 ^ Rc- To find r'g in this limit, one should 
solve the following auxiliary problem. Collect all closed 
equipotential contours of size of order ^5 S> (i within the 
area ^5X^5- They form a "bundle" of width 



ws ~ d{d/vdr'sf/' ^ d{d/is?/^ 



(31) 



The characteristic area covered by this bundle is S{ws) ~ 
Lsws ^ C|. Now enlarge the area by adding all points 
that are within a distance A > ws of the initial bun- 
dle. Doing so we get a new bundle that occupies an 
area ^(A). The question is what is S'(A) at A '--^ i?c for 
> > d. Clearly, S'(A) ~ LsA as long as A < d. 
However, at A ^ d, the area should exhibit a scaling 
behavior S'(A) ~ Lsd{A/d)^ with a nontrivial exponent 
X. To find X, notice that at A ^ ^5 we should have 
'S'(Cs) ~ S.S- Using the equation Ls ~ d{£,s/dy^'^ we 
thus get X — l/A. It follows that the area S{Rc) scales 
at i?c ^ d as Lsd{Rc/dy^^ [which can be represented as 
S,siBc/£,sy^^ to see that most of space within the area 
£,S X £,S is left empty]. Since (5 <C a, it is clear that if 
there is an AD in this area, it will be inevitably hit by 
the cyclotron orbit. Therefore, the length Ls obeys the 
equation 



nSiRc) - 1 , 
which yields Ls ^ {nRcy^{Rc/df/^ and 



TS 



3/4 



(32) 



(33) 



We see that Eqs. (|30|),(33) match each other at Rc ^ 
^ ^ On the other hand, Eqs. (|^),(|3l) match at 
Rc - d. 

The above derivation oi r'g a.t 5 /a ^ 1 shows that the 
increase of the scattering time as compared to the "Drude 
time" Ts is due to multiple passages of the cyclotron orbit 
through the area ax a. Clearly, at Rc/d <C 1 no renormal- 
ization of the scattering time occurs as long as 5/a S> 1. 
However, for Rc/d 3> 1 the scattering time is renormal- 
ized with increasing B already at some 5/a^ 1. Indeed, 
let 5/a 3> 1 and consider the case Ls ^ d, <C Rc- The 
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drifting cyclotron orbit experiences Ls/^s ^ 1 returns 
to the area d x d and each time it probes the fraction 
a/ 6 <C 1 of space within this area. One sees that if the 
product of the two factors {Ls/^s){0'/S) » 1, then the 
colhsion time is much larger than the Drude time ts and 
obeys Eq. (|9|), which yields Eq. (|3^). Therefore, it is 



only when 6/a '3> (nRcd) '^1^ ^ 1 that r'g = ts- Simi- 
larly, if > Rc, Eq.^) is vahd for aU 6 /a < {Rc/df^. 

Inspection of Eqs. (|27|),(p0|),([3^) shows that the depen- 
dence of Tg/rs on B is parametrized by two parameters. 
In addition to p [Eq. ([l^)], it is convenient to introduce 
the parameter 



rj — nd^p , 



(34) 



which can be rewritten as Vod/ej^a, where Vq is a charac- 
teristic amplitude of fluctuations of the long-range poten- 
tial, ep the Fermi energy. The meaning of this parameter 



is that it describes the position of the crossover S/a ^ 1 
in terms of the ratio d/Rc- Specifically, if r/ <C 1, the 



crossover occurs at d/ Rc 



d/Rc 



1/2 



2/3 



whereas if 77 3> 1, at 



Note that the crossover point at 3> 1 



corresponds to Rc ^ a [namely Rc/a ~ (dZi/a^)^/^], 
however large rj is. We thus have 



(35) 



By changing rj with p held constant, we change the con- 



centration n at a fixed mean free path Is- In Sec. IIIB, 
we have already studied the limit 77 — > cxd: in that case 
the scattering time is not renormalized, so that r'g /ts = 1 
for aU B independently of p. However, Eqs. (|^) , (|ol) , (H) 
tell us that at any finite 77 there exists a magnetic field 
above which t^/ts starts to grow with increasing B. The 
function h{x,p,r]) which describes this growth reads 



h{x,p,ri) 



1 — 1 
a; 77 

j;3/2y^-l 

2,9/4^3/4^-7/4 

2;3/4j^-l 



» max{77^/^, 1} ; 



,2/3 



< X < min{7775"\ 1} , < 77 < 1 ; 
max{77p~^, if I'^p-'^l'^^ ^ X <C rj^l'^p^^l'^ , 77 ^ min{p,p^'^/^} 
max{77'*/^, ri^l'^p-'^l'^^ < a; < 1 , 77 < min{p, 1} . 



(36) 



The behavior of t^/ts given by Eqs. (|35|),(|36D is illus- 
trated in Fig. 1^. One sees that in the limit of large B 
the collision time grows as i?^, whatever the parameters 
p and 77. If 77 ^ 1, which corresponds to a sufficiently 
large concentration of ADs, this B^ growth matches the 
Drude result directly [curves labeled by (i) in Figs. |^a,b]. 
In a more dilute array of ADs, there appear intermediate 
regimes, which proliferate as 77 is decreased [curves (ii)- 
(iv)]. Note that the smaller 77 for a given p, the sooner 
the renormalization of r'g starts with increasing B. 

In Fig. we marked the point d/Rc ^ {d/lhY^^ 
which corresponds to the crossover between the diffu- 
sion and drift regimes [smaller and larger d/Rc respec- 
tively, see Eq. (p^)]. Equations (|35|),(^ describe the 
drift regime. However, as illustrated in Fig. |^ by the 
uppermost (dashed) curve, if the concentration n is suf- 
ficiently low, the collision time is renormalized already 
in the diffusive regime. Indeed, for small 77,, t^/tq calcu- 
lated for drifting electrons will be large at the crossover 
point to the diffusive regime. A similar effect takes place 
for large p as well (not shown in Fig. Qb). In the dif- 
fusive regime, Tg/rs is still given by Eq. (|29|), the only 
difference is that now in place of Eq. ( p8|) one should take 



^s-S{vft's/Rc)'/^ , 



(37) 



which describes the diffusive motion of the cyclotron or- 
bit. Substituting this expression for 1^5 in Eq. (E9h yields 



II 1 

Is {nRlf 



(38) 



r 



i.e., the collision time starts to grow as i?"* before enter- 
ing the drift regime, as shown in Fig. |7^ by the dashed 
line. In the diffusive regime, Pxx{B) is related to r'g by 

Pxx 

{B)/po Ts/Tg, which gives 



Pxx{B) Is 

Tj 



Po 



{nRl 



(39) 



As follows from Eq. (|3^) , the crossover to the B'^ behav- 
ior with increasing B occurs at 



d/Rc^ind'Y'^ils/lLY'^ . 



(40) 



Comparing Eq. ( ^ ) with (d/l^)^/^ we conclude that the 
crossover to the drift regime takes place at larger B if 
nd'^ ^ {d/lLY/^[lL/lsY^'^ ■ If, however, the concentra- 
tion n is high enough, so that the opposite inequality 
is met, the region of validity of Eqs. (|3^),(p9|) shrinks 
away and no renormalization of ts occurs in the diffusive 
regime. 

Note a sharp change in the behavior of Tg/rs at 
the crossover between the drift and diffusion regimes 
(Fig. I^a). The mismatch is due to the difference in 
the fractal dimensionality of extended trajectories in the 
two regimes: self-avoiding drift trajectories percolate in a 
super-diffusive way and therefore explore the area faster. 
The sharp crossover has the form of an exponential falloff 
of Tg/rg, governed by rapidly developing adiabaticity of 
electron motion. If one compares the times r'g obtained 
for the two (diffusion and drift) regimes close to the 
crossover point d/Rc ^ {d/lL)^^^, one can see that for 



10 



small n, namely for ntP <C {(I/Il)^/'' {Il/IsY^"^ , the ratio 
Tg/rs S> 1 on both sides. The dashed line in Fig. |^ 
illustrates just this case. It is possible, however, that 
n falls into the intermediate range {d/lLf^'' (Il/Is)'^^'' ^ 
n(P <C (d/^L)^^^(^L/^s)^^^, in which case Tg changes with 
increasing B in the following way: it first starts to grow 
in the diffusive regime, then, after the crossover into the 
drift regime, returns to the unrenormalized value ts, and 
only with further increasing B begins to grow again. 



Ts/Ts 



a) 





9/4/ 9/4/ 











(i) 



3/2/ 





1 


d/Rc 






/(iii) 


b) 




2 / 

/ / (ii) 
V 

/ / (i) 


9/4/ 


3/4^^ 


V 









1 



d/Rc 



FIG. 7. Schematic behavior of r'g/rs as a function of d/Rc 
on a log- log scale for a) p <C 1 and b) p ^ 1. The numbers 
denote the exponent of the power-law dependence of rg /ts on 
B. Different curves in Fig. ^ illustrate how the dependence 
of t's/ts on B for a fixed p is modified in different ranges 
of r]-. (i) 77 > 1; (ii) p < t? < 1; (iii) -q p; (iv) 

77 < Similarly in Fig. |b: (i) 77 > 1; (ii) p"^/^ < 77 < 1; 
(iii) T) <^ p~^^^ . The dashed line in Fig. Ma shows the behav- 
ior of t's/ts for r] so small that the renormahzation of t's/ts 
starts with increasing B already in the diffusive regime (see 
the text). 

At this point it is worth emphasizing once more that, 
while in the diffusion regime all electrons behave in a 
similar way and given by Eq. (|3^) is characteristic 
to all electrons, upon crossover into the drift regime the 
electrons find themselves divided into different groups. 



characterized by different collision times (the groups 
are mixed up only due to slow nonadiabatic dynamics). 
Specifically, in Eq. ( pO[ ) is the collision time for elec- 
trons that move along extended (percolative) drift tra- 
jectories. This time is renormalized even for S ^ a. On 
the other hand, electrons that upon crossover to the drift 
regime find themselves on typical drift trajectories of size 
~ d either do not collide with ADs at all, and for them the 
collision time is infinite (actually, provided they experi- 
ence nonadiabatic transitions, the collision time is finite 
but exponentially large, as shown in Sec. 



HIE) 



characterized by the unrenor maliz ed collision time rg (at 
(5 ^ a), as explained in Sec. Ill C. Hence, the renormal- 



ization of the collision time that we have analyzed in this 
section will affect the percol ative c ontribution to the MR, 
which we will study in Sec. IVD. For 5 <C a, the renor- 
mahzation of Ts will strongly affect also the short-scale 
dynamics of electrons residing on typical trajectories of 



size d, which will be considered in Sec. [V E 



Prr(B) 




FIG. 8. Schematic behavior of pxx{B) as a function of 



,1/3 



Dif- 



d/Rc in the diffusive regime d/Rc <^ {d/lL 
ferent curves illustrate the dependence of pxx{B) on B 

(i) nd^ » (d//t)2/3(/L/;s)^/2. 

(iii) 



in different ranges of n: 
(ii) {d/hf^^ < nd^ 
nd^ « (d/hf^^ 



« {d/lLf^'{lL/ls) 



Now let us return to the diffusive regime. Clearly, 
Eq. (^) is only valid as long as vft'^ <C II, which is 
rewritten as nR^ ^1. If the opposite limit, nR^ <C 1, 
is realized with increasing B still in the diffusive regime 
[which is the case for a very dilute array of ADs, namely 
for nd^ <C (d/Zi)^/'^], in this limit the scattering on ADs 
stops playing any role for diffusive electrons. The colli- 
sion time is then given by I l /vp and 



Pxx{B)/po = Is/h 



(41) 



does not depend on _B, whereas t'^/ts keeps growing with 
increasing B. Yet, the scattering on ADs will become rel- 
evant again with further increasing B, once the system 
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crosses over into the drift regime (where the scattering 
on ADs wiU prevent the adiabatic locahzation from com- 
pletely suppressing the MR, as explained above). This 
behavior corresponds to the case when the crossover to 
the diffusive regime with decreasing B occurs not in the 
region of i?^/"* behavior, as shown in Fig. but in the 
^3/4 j.ggiQ]-^^ xhe overall behavior of the MR in the dif- 
fusive regime is illustrated in Fig. ^. We consider the 
effect of the renormalization of rg for diffusive electrons 
in detail elsewhere [3^ ]. 

As mentioned above, another effect of small 6 /a is a 
renormalization of the hopping length for the diffusive 
dynamics across the drift trajectory due to the scatter- 
ing on ADs, which we will consider in more detail below. 
Since both the effective scattering time and the hopping 
length are now modified, the shape of the MR is no longer 
parametrized by a single parameter. Specifically, we can 
replace / in Eq. (^^ by a new function g which is given 
by 



gix,p) r/ h^' P 



TS , 



(42) 



in the "one-hop" regime [cf. Eqs. ( [l5|) , (|l6|)] , while in the 
advection-diffusion regime [cf. Eq. (|l7|)] we now have 



( \ f Rh 



6/13 



fix 



P— 



(43) 



A crossover between the two regimes occurs at vaTg ~ 
L{Rfi), where L{w) is defined in Eq. (Q). The ratio t'^/ts 
is given by Eqs. (p5|),(|3^). The ratio RyJ will be cal- 



culated in Sees. IVB IVC 



B. Skipping orbits 



In Sec. [VA, we derived general expressions for the 



percolative MR for small (5/a in terms of r'g and Rh and 
calculated the effective scattering time r'g. Let us now 
evaluate Rh- The scattering problem for a cyclotron or- 
bit that collides with a hard disc is nontrivial dX 5/a <^ 1. 
To begin with, notice that at 5/a <C 1 the drifting cy- 
clotron orbit first hits the disc boundary at a small angle 
9i <^ 1 (see Fig. ^. A simple geometric consideration 
yields a characteristic 9i for a particle incident on the 
disc with a drift-shift vector 5 (i.e., with a drift velocity 
8vf/2ttRc): 



1/2 



Rc + a 
Rc 



1/2 



(44) 



where e is the unit vector normal to the surface of the 
disc at the point of the collision. One sees that the an- 
gle of incidence vanishes when we let ^ ^ 0, so that 
the skipping cyclotron orbit in effect starts to "roll over" 
the disc. Clearly, if it were not for the drift during the 



rollover, the collision angle 9 would be conserved, being 
the same each time the particle returns to the disc after 
one cyclotron revolution. Therefore, although the drift 
is slow, in the sense that the typical hopping length for 
the guiding center after one collision ~ Rc9 » 6, it is 
because of the drift during the skipping process that the 
particle eventually breaks away from the AD. The length 
Rh is then understood as a shift of the guiding center at 
the point of the breakaway with respect to the equipo- 
tential contour along which it was drifting right before 
the first hit. 




FIG. 9. Geometry of scattering of a skipping cyclotron or- 
bit on an antidot: 6 is the angle of incidence, (j> the polar angle 
of the point at which the collision occurs, ^j) the incremental 
increase of the angle between two consecutive collisions. 

Since 5/a <^ 1, we can treat the drift during the 
rollover perturbatively. The system of equations that 
describe the skipping in the absence of drift is given by 



K+i = 4>7i + i>{On) ; 

9n+l = 9n , 



(45) 
(46) 



where 0„ is the polar angle (along the surface of the disc) 
defining the point of the n-th collision, 0„ the collision 
angle, see Fig. ^. The function ?/'(6') obeys the equation 
i?c sin(6' — ^/2) = asm{ip/2), which reduces to the linear 
relation 



2Rc 



Rr 



(47) 



in the limit of small 9. Equation (^) says 9 is the integral 
of motion for the skipping process without drift. In the 
presence of drift, 0„ acquires an n-dcpendcnt correction 
A9n = 9n+i — 9n. To first order in S the correction reads 



A9r^~- 



Rc + a Sne{4>n) 

RnU 9n 



(48) 



where 5n is the drift shift between the n-th and (n-l-l)-st 
collisions, e(0„) the unit vector perpendicular to the sur- 
face of the disc at the point of the n-th collision. Trans- 
forming to the continuous limit we get the differential 
equation 
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2 R^ + a 



d8^ 



(49) 



whose solution, after substituting ip{9) from Eq. (47) 
yields 



^(0) = 0^(^.)-7r 



3 [Rc + a 



2a V Rc 



d(t)'8{(t)')e{(t>') . (50) 



The function 5(0) in Eqs. (g|),(|§ gives the drift shift 
for the cyclotron orbit whose guiding center is a distance 
Rc+a from the center of the disc in the direction specified 
by the angle 4>. In fact, the integration in Eq. (^^ runs 
along the guiding-center trajectory during the rollover, 
which is the contour p = R^. + a, where p is the radius 
vector counted from the center of the disc. Note that the 
integral term in Eq. (|o|) is bounded from above by ^ S/a 
at Rc ^ a and by 6a/ R^ ^ a/ {dli,)^/'^ otherwise. In both 
limits the maximum value of <IC 1, which justifies the 
linearization ( ^7| ) and our using of the term "rollover" . 

If we take (pi in Eq. ( |50| ) equal to the polar angle 
at which the cyclotron orbit hits the disc for the first 
time, then 0{<j)i) should be put to zero in the continu- 
ous approximation. With the same accuracy the angle 
4>f at which the breakaway occurs satisfies the condition 
9{(j)f) — 0. According to Eq. (pO|), we can recast the 
latter condition as 



dp <5(0)e(0) = 



(51) 



This equation yields cf>f as a function of pi and, conse- 
quently, enables us to determine the shift Rh- Let us 
introduce an effective random potential V{p) as the av- 
erage of the real potential over the cyclotron orbit with 
the guiding center at the point p [at Rc <^ d the two po- 
tentials almost coincide, but at Rc':^ d a typical ampli- 
tude of fluctuations of the effective potential V{p) with 
the same correlation radius d is obviously ~ {Rc/dY^^ 
times smaller] . The drift occurs along equipotential lines 
of V{p). Since the shift 6{p) oc W{p) x e^, where 
is the unit vector along the magnetic field, Eq. ( |5l| ) is 
rewritten as 



4'j 



{VV X e,) ■ e{p) = 



(52) 



which, for the integration along the arc p = Rc+a, finally 
gives 



V{p,Pf)-V{p,cl),)=0 



(53) 



We thus see that in the limit of small 6/a the cyclotron 
orbit, having skipped along the surface of the hard disc, 
breaks away on the equipotential contour with the same 
V as it had before hitting the disc, i.e.. 



in the continuous approximation. Note that for the drift 
in a homogeneous electric field, when 5{p) = const, this 
result follows straightforwardly from symmetry of the 
scattering problem. What Eq. ( ^3| ) tells us is that, re- 
markably, Rh vanishes in the case of varying V{p) as 
well. 

Let us compare the above picture with a familiar ex- 
ample of adiabaticity of scattering on a smooth inhomo- 
gcneity V{p): in that case, the vanishing of Rh simply 
means that the guiding center drifts along a locally per- 
turbed equipotential line of V{p). Naively, one might 
think that a collision with the disc destroys the adia- 
baticity since the impenetrable hard disc makes a part 
of the equipotential line of V{p) inaccessible. However, 
as follows from Eq. (53), the skipping of the cyclotron 
orbit around the disc goes on adiabatically, provided S/a 
is infinitesimally small, and the result Rh — still holds. 



C. Nonadiabatic skipping 



Rh = 



(54) 



The zero result for Rh was obtained in Sec. [VB by 
treating the drift during the rollover perturbatively, to 
first order in S/a, and by taking the continuous limit. To 
find Rh, we should now relax this approximation. Before 
doing so, it is worthwhile to recall how Rh behaves for 
scattering on a smooth inhomogeneity. In that problem, 
it is known that taking higher gradient terms into account 
does not lead to a finite Rh and, in fact, the problem 
of finding Rh does not allow for any perturbative solu- 
tion that could be expanded in powers of the parameter 
S/d <ti 1, where d is a characteristic size of the inhomo- 
geneity. Specifically, for a smooth inhomogeneity, Rh is 
exponentially small at i5/d <C 1 (see, e.g., [Q for a solu- 
tion of the scattering problem and references therein) . As 
we have shown above, the case of a hard disc placed in a 
smoothly varying environment is similar in that the scat- 
tering is also almost adiabatic at (5/a <C 1. A question 
then arises if the nonadiabatic scattering that leads to a 
finite Rh is also exponentially suppressed. The answer 
is no, since there is an important difference between the 
two cases. Namely, the approximation within which the 
skipping of the cyclotron orbit can be considered as a con- 
tinuous adiabatic process [Eq. (49)] fails completely near 
the points p = pi and p = pf. Indeed, near these points 
dynamics of the collision angle 9(p) is nonadiabatic since 
9 is close to zero, so that A6'„ is of order 9n itself and 
the expansion (^) is not valid any more. A finite shift 
Rh is therefore due to the discreteness and incommensu- 
rability of the skipping along the sharp boundary of the 
hard disc. It is given by the elementary (associated with 
a single collision) hopping length of the guiding center 
near p = (t>i,f^ which is Rhi ^ {Rc + a)'ip{9i). Substi- 
tuting 9i from Eq. (^^ we finally get for a characteristic 
amplitude of the shift 
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Rhi ^ Rc 



S Rc 



Rc 



1/2 



(55) 



We see that, in contrast to the case of a smooth inhomo- 
geneity, R^i scales as a power of S, namely R^i oc S^^^. 
Note also that the characteristic scale d, on which the 
smooth V{p) changes, appears in Eq. (js^) only through 
the drift shift 5. 

Equation (|55| ) describes a single scattering in which 
the particle breaks away from the disc along an equipo- 
tential with V almost equal to that of the equipotential 
along which it is incident on the disc. Using Eq. (55), 
dynamics of skipping orbits in the limit Rc <^ d can be 
understood quite straightforwardly. Since in this limit 
the drift is almost homogeneous in the course of skip- 
ping, for a given V there is typically only one equipoten- 
tial line that crosses the guiding-center trajectory during 
the rollover. Accordingly, if it were not for the small 
shift (|5^), the particle would simply continue to drift af- 
ter the rollover along the same equipotential line. The 
picture becomes far more complicated at large Rc ^ d 
since in that case there are many equipotential lines with 
the same V that intersect the guiding-center trajectory 
corresponding to the rollover. 




FIG. 10. Schematic picture of scattering of a cyclotron or- 
bit on an antidot a,t S <ti a, Rc ^ d and R^s ^ ws- The 
thick solid line with arrows is the guiding-center trajectory. 
The dashed lines are equipotential contours which cannot be 
accessed by the guiding center. The small shaded circle in the 
center of the figure shows the antidot. The thin circle of radius 
Rc + a around it is the boundary of the area impenetrable for 
the guiding center. The parts of the guiding-center trajectory 
that coincide with this boundary correspond to skipping of 
the cyclotron orbit, which alternate with parts corresponding 
to drift. The drift occurs between each breakaway from the 
antidot and consecutive return to it. 

Let Rc ^ d. In this limit, one should distinguish two 
regimes according to whether the typical hopping length 
after one coUision Rhi ~ RdS/aY^^ [Eq. (|55|)] is larger 



or smaller than d. Consider first the case Rhi <C d, i.e., 
let d <C i?c d{a/6y^^. The characteristic number of 
equipotential lines with the same V that cross the circle 
of radius p = Rc + a around the disc (which is the guid- 
ing center trajectory in the course of skipping) is of order 
Rc/d ^ 1. Clearly, the direction of drift (to or away from 
the surface of the disc) alternates during the skipping. It 
follows that the particle which started skipping will break 
away along the equipotential line that is the first to cross 
the guiding-center trajectory corresponding to the skip- 
ping. Yet, since the equipotential lines are closed loops 
(typically of size d), the particle will come full circle and 
return to the disc (see Fig. ^o|). Then the process will re- 
peat itself with other equipotential lines along the surface 
of the disc. We have assumed, however, that the particle 
is incident on the disc along a percolating extended tra- 
jectory. Therefore, the multiple returns will stop when 
the particle picks up this trajectory. 

Let us evaluate the total shift R^^ with which the par- 
ticle will finally break away. Elementary shifts for the 
repeating collisions are uncorrelated with each other, so 
that 



d2 d2 p 

J^c 



dKKW{K) 



(56) 



where iy(A) is the probability density for the drift tra- 
jectory that broke away from the disc to hit it again, i.e., 
to cross the arc p = Rc + a, for the first time after the 
breakaway at a distance A from the starting point. It is 
instructive to map the problem of finding VF(A) onto a 
more conventional one by noting that the power-law scal- 
ing of W{K) describes how a deposition rate for particles 
emitted by a point source and moving in two dimensions 
in the presence of an absorbing line falls off with increas- 
ing distance A along this line. If the particles would 
experience an uncorrelated diffusion with the elementary 
step ^ d, then it is straightforward to see, by solving 
the diffusion equation with the absorbing boundary, that 
iy(A) ~ d/A^. In fact, one can show, by introducing a 
scale-dependent diffusion coefficient which describes the 
drift, that this result holds for the drifting particles as 
well, i.e., W{K) and P(^) [Eq. (|l9|)] have the same scahng 
behavior. Hence, the integral in Eq. ( |56| ) logarithmically 
diverges and the typical number of collisions before the 
final breakaway is Rc/d\n{Rc/d). It follows that 



R 



hT. 



Rc 



Rc 



-,1/2 



a d\n{Rc/d) 



(57) 



Clearly, the necessary condition for the above derivation 
of RhT, to be valid is Rh^ <^ d <^ Rc, which means 
d<t: Rc<^ d{a/S)^^^ ln^/^(a/(5). It is worthwhile to men- 
tion that if Rc ^ d{a/Sy^^ the characteristic hopping 
length for the skipping cyclotron orbit Rctpmax ^ d. 
Here ipmax (<5/a)^/"^ is the maximum scattering angle 
for a single-run skipping at i?c ^ a [see Eqs. (E7|),(|5r 
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Now let Rhi be still smaller than d but let Rhu ^ d, 
i.e., consider the interval d{a/S)'^^^lii^^^{a/S) < i?c < 
d{a/Sy^^. An essential difference appears in this regime: 
the shift accumulated through multiple breakaways and 
returns exceeds d before the rollover is finished. At this 
point it is worth recalling that the elementary shifts for 
each breakaway are accompanied by changes of V corre- 
sponding to drift trajectories. Accordingly, when the ac- 
cumulated shift gets larger than d, the initial and current 
values of V become uncorrelated with each other. It is ev- 
ident that we should treat this case separately, since now 
one cannot identify Rhu given by Eq. ( p7[ ) with a shift 
with which the particle has finally broken away never to 
return. This conclusion becomes even clearer when not 
only RhT, but also Rhi [Eq. (|55|)] is larger than d, i.e., 
when Rc S> d{a/5y^^. In the latter case, the memory 
about the initial value of V is lost already after one col- 
lision. 




FIG. 11. Schematic picture of scattering of a cyclotron or- 
bit on antidots a.t S <^ a, Rc <^ d and Rhi ^ ws- The line 
with arrows shows the guiding-center trajectory. The thicker 
line is the trajectory that takes the particle from one antidot 
to another (two of which are shown by shaded circles) . Until 
the guiding center picks up this trajectory the particle is stuck 
to an antidot and keeps colliding with it, by making long drift 
excursions between the collisions. For Rc ^ d, the scattering 
process can be visualized as a "combination" of Figs. p^pd] . 

The picture emerging in the limit RhYi 3> d signals that 
one should also fine-tune the derivation of Eqs. (^),(57). 
Namely, these equations should be supplemented with 
the condition under which the problem of finding Rh can 
in fact be formulated as the scattering problem for a sin- 
gle disc. Indeed, for a single disc, one actually cannot 
specify with which accuracy the particle should hit the 
vicinity of the percolating trajectory so as to be able 
to break away from the disc. To answer this question, 
we need to consider a scattering problem for two discs. 
Specifically, consider two discs separated by a typical dis- 
tance Ls ~ v^t'c; measured along the percolating trajec- 
tory (Fig. |l^). The two discs are connected by a bun- 
dle of drift trajectories of width ws [Eq. (^l|)]. One 
sees that the particle will break away from one disc and 
get through to the other after a single rollover only if 
max{i?;ii, i?/i5]} is within this width. If this condition is 
not satisfied, the cyclotron orbit will keep going round 
the disc until it happens that its guiding center hits the 



strip of width ws at the same time when vanishes. It 
follows that at Rc d the shift is given by Eq. ( [55| ) 
only if Rhi <C ws- If, however, R^i ^ ^5, it takes typ- 
ically Nj. ~ Rfii/ws revolutions around the disc before 
the particle breaks away and 



Rh ^ Ws 



(58) 



in this limit. Similarly, ed, Rc^ d the shift R^ is given by 
Eq. (|5^) only if Rhs «C ws; otherwise, Rh ~ ws (the pic- 
ture can be visualized as a "combination" of Figs. |lo| , pT| ). 
For Rc ^ d, the number of revolutions around the disc 
Nr, necessary for the simultaneous tuning of V and 9 cor- 
responding to the breakaway, is different in three different 
regimes. Namely, iV^ ^ RhT./ws for ws ^ RhT, ^ d, in- 
dependently of the ratio Rhi/ws- With increasing Rh^, 
if Rhs ^ d, we have Nr ^ d/ws as long as Rhi <C d, 
and Nr ^ Rhi/ws otherwise. In particular, one sees that 
if Rhi: S> d, it always takes many revolutions before the 
breakaway occurs and the shift Rh is always determined 
by ws- 



D. Percolation in a dilute antidot array 



In Sec. IV C , we discussed scattering of a drifting cy- 
clotron orbit on a single AD at (5/a <C 1. Consider now 
the MR at (5/a ^ 1. Let us start with the limit of a sm all 
concentration of ADs n ^ 0. As outlined in Sec. [IB, in 



the absence of ADs the conductivity at large B is only 
due to the exponentially weak nonadiabatic scattering on 
long-range disorder. Let us neglect this additional con- 
tribution to Pxx{B) and calculate the contribution that 
is due to the scattering on very rare hard discs. We thus 
seek a term in Pxx{B) which is proportional to a power 
of n at small n. Let Rc <^ d. Since at n — > the bun- 
dle of drift trajectories that connect two discs becomes 
infinitcsimally narrow [ws oc n'^/'' at n — > according to 
Eqs. (p7|) , (pT])] , the particle sticks to a disc for a long time 
Tst until it picks up a trajectory that is extended enough 
to take it to another disc. The number of the unsuc- 
cessful attempts to break away from the disc is given by 
Nr ~ Rhi/ws 3> 1. To evaluate r^t, notice that the stick- 
ing time is determined by a slow drift along the closed 
loops that repeatedly return the particle back to the disc, 
not by the fast skipping in between. However, most of 
the attempts end up in quick returns to the vicinity of 
the point of the first collision and so give only a small 
contribution to Tst ■ We estimate Tst as 



Tst 



^" dw L{w) 

Ws Vd 



(59) 



where L{w) is given by Eq. (|^). The integral is deter- 
mined by w Ws, which yields Tst Tg. We thus see 
that the effective scattering time between collisions with 
different discs is given by the drift time between the discs. 
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Accordingly, the total distance passed by the particle in 
the multiple rollovers is typically of order Ls- In fact, 
this conclusion holds for the case Rc ^ d a.s well. The 
only difference is that at large Rc the particle experiences 
multiple breakaways and returns on the scale of a single 
rollover. However, Tst does not renormalize the charac- 
teristic time between collisions with different discs, which 
is given by r'g, similarly to the case Rc ^ d. 

We are now in a position to calculate the AD-induccd 
contribution to the MR in the limit of small n. Accord- 
ing to the picture above, the scattering on a given AD is 
over when the particle hits another AD separated by a 
distance ^ dlvdT'g/d)^^"^ , and the characteristic time 
between two scatterings is r'g. As long as we neglect the 
nonadiabatic corrections to the drift in the long-range 
potential, there is a clear separation between localized 
and extended drift trajectories. Only a small fraction of 
trajectories get delocalized by means of collisions with 
ADs, namely ^ wsLs/^.g <C 1. The diffusion coefficient 
of particles residing on these trajectories is ~ Vd(^/Ls- 
We thus see that the macroscopic diffusion coefficient is 
given by Eq. (|^) with a rescaled scattering time ts — > t\ 
which yields the MR obeying Eq. (^. 

It is worth noting that the dynamics of particles is es- 

(Sec 



sentially different at 



oo 



[II B) 



and n — > 0, 

despite being described by similar equations. In the for- 
mer case, there is a strong exchange, caused by collisions 
with ADs and governed by a detailed balance of scatter- 
ing processes, between the stream of fast particles which 
follow the links of the percolation network and the reser- 
voir of "quasilocalized" particles which stick for a long 
time to within the critical cells of the network. By con- 
trast, at n — !■ there is no such exchange and the drift 
trajectories within the cells of the percolation network 
are strictly localized. However, in both cases the MR is 
determined by the fast particles moving along the links of 
the percolation network, which is why it does not matter 
if the particles moving inside the critical cells are local- 
ized or not. Put another way, although the total number 
of delocalized particles decreases at n ^ 0, this effect is 
compensated by more frequent crossings, due to collisions 
with ADs, of the percolative drift trajectory by particles 
that remain delocalized. 

Substituting Eqs. (|3),(|3|) into Eq. (El we get 



for Rc^ d and 



Pxx{B) _ / d 

'^^ ' \Rc 



P^^jB) ^ g4/7 3/7 

Po \ Rc 



39/28 



4/7 



(60) 



(61) 



in the opposite limit. According to Eqs. (|60[) , (p3l|) , pxx oc 
n'^/'' at n 0. 

One sees from Eq. (^) that pxx taken in the limit 
n — > behaves at large B as B^/"^ . As compared to 



Eq. ( p7| ) , which describes the asymptotic behavior in the 
hydrodynamic regime n — s- oo and gives pxx B^^/^^, 
the divergence of Pxx with increasing B is weakened, but 
is still present. However, neither of the two limiting cases 
( p7| ) , (|6^) describes the asymptotics of the MR for i? — + cx3 
at a given finite n, which we discuss below. 

Let us turn to the percolative MR in a denser array 
of ADs. Increasing n yields wider links of the percola- 
tion network, so that eventually wg becomes larger than 
the elementary shift ma.yi{Rhi-,RhT} (Sec. [V C| ). This 
transport regime is described by the advection-diffusion 
equation (17) modified according to Eq. (142 
^,(57) for Rhi 



Eqs. (^ 
Tg in Eq. (43) gives 



and and Eqs. 




, plO/13^6/13 

Po V Rc 



1/52 



ln-3/13 



Rc 



for Rc^ d and 

Pxx{B 



PO 



R \ 2/13 / p I „ \ 3/13 
pio/13^6/13 I / nc + a 



(62) 



(63) 



for Rc < d. 

The condition at which the "one-hop" percolation 
[Eqs. (|60|) , (|6l|)] crosses over with increasing B into the 
advection-diffusion regime [Eqs. (|6^),(|63|)] is given by 
somewhat cumbersome formulas: the crossover occurs 
at d/Rc - (/r/)3/i25ln"2^/^2^(p6r;) for p^r/ < 1, at 
d/Rc - (/ry)i/22 for 1 < (^77)1/22 < d/a, and at 
d/Rc - (/r;)i/i5(a/d)7/i5 for (^77)1/22 > d/a. 

Equation ( |63| ) tells us that the asymptotics of the MR 
at B — > 00 is Pxx oc iJi/i'^, i.e., the MR diverges as a 
power law in the limit of large B. However, this diver- 
gence is so weak that from a practical point of view it is 
indistinguishable from a saturation of pxx (B) . Nonethe- 
less, it is a remarkable fact that in the extreme i? — > 00 
the MR in the presence of both ADs and a long-range 
potential does not go to zero, in contrast to Eqs. (|l|),(|j), 
which predict vanishing of pxx when only one type of 
disorder is present. 



E. Rosette orbits in the presence of weak drift 



In Sec. IV D, we considered the percolative contribu- 
tion to the MR at small 6 /a. Now we proceed to the 
"short-scale" contribution associated with rosette states 



(Sees. [IA,[IIC). Let us analyze how the Lorcntz-gas be- 
havior is restored with decreasing strength of the smooth 
disorder. Since we deal with the case a/d ^ 1, the 
Lorentz-model limit is achieved for a very weak long- 
range potential, such that at Rc/ls ^ Ij when the 
falloff (|l|) starts, the scattering on the long-range po- 
tential is already strongly adiabatic. The condition of 
adiabaticity at Rc/ls ~ 1 and Is/d'^ 1 translates into 
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Ih/d » {Is/d)^. As we will see below, the condition 
of the Lorentz-gas falloff Px{B) oc not being domi- 
nated by the contribution of "hopping rings" (Sec. Ill C ) 
is much stronger, namely 



lLld-> {lsldf{d/aY 



(64) 



The solution of the scattering problem for a single disc 
in Sec. IV B| shows that at (5/a ^ 1 there is a well-defined 
separatrix Omax ^ 1 for the angle of incidence 6 in the 
phase space (9,(1)), which divides skipping cyclotron or- 
bits into two groups, delocalized and localized. Namely, 
trajectories that will finally break away from the disc 
belong to the part of the phase space with 9 < 9max, 
whereas the region 9 > 9„iax is filled with those that will 
never escape. In other words, almost the whole phase 
space is filled with bound states. According to Eq. (50), 
the critical angle 



{d/ay/^[iR, + a)/R, 



12/3 



(65) 



Indeed, a drifting cyclotron orbit that is incident on the 
disc both hits it for the first time and finally breaks away 

9 < 6i [Eq. (^)], while reaching a maximum 9, which 
depends on the incident parameter but does not exceed 
9max given by Eq. (|6^), during a rollover in between. 

Let us show that a particle which starts with some 
9 > 9max will not be able to break away. What is im- 
portant to us is that in order to break away the par- 
ticle has to decrease 9 down to 9 ^ 9i. Suppose this 
might happen and the particle has escaped. Then we 
could consider a complementary scattering problem by 
reversing time and sending the particle that has broken 



away back to the disc. However, we know from Sec. IV B 



that for this scattering problem 9 will never exceed 9max- 
Hence, we come to a contradiction which shows that the 
initial assumption about the possibility of a breakaway 
cannot be realized. We thus see that trajectories that 
started to skip along the disc with 9 > 9max can never 
cross the boundary 9 = 9max and so will remain bound 
to the disc. The fact that there exist bound states for 
magnetized electrons interacting with a single hard disc 
in the case of a homogeneous in-plane electric field was 
observed in the numerical simulation [^7| . 

The existence of the separatrix 9max means that the 
contribution to Pxx{B) of rosette states with 9 > 9max is 
given by 



Pxx{B)/po ^ Rc/h 



(66) 



for aU 5/a < 1 [cf. Eq. 

Let us now compare two terms in Pxx{B) that are as- 
sociated with the short-scale diffusion at finite 5 /a. One 
described by Eq. ( p2|), i s due to the hopping of "rings'' 
introduced in Sec. IIIC. The other, given by Eq. (66) 



is due to the hopping of rosette states. Note a similarity 
between the two mechanisms of transport: in both cases 



a particle interacts with a disc many times before chang- 
ing to another disc. The comparison shows that the con- 
centration of particles participating in the hopping-ring 
transport is much larger than that of rosette states. It is 
clear, however, that if we send ^ — > 0, the hopping rings 
should not contribute to pxx{B), which will be given by 
Eq. (^. It follows that there should exist yet another, 
intermediate regime of hopping, associated with the evo- 
lution of the hopping-ring transport with decreasing S/a. 
To describe the latter, notice first of all that Eq. ( |2^ ) 
stops to be valid already at some 6 /a » 1. Indeed, the 
logarithmic factor in Eq. (|2^) is associated with the drift 
along closed loops of size ^ < (ni?c)~^- Therefore, the 
derivation of Eq. (^2|) in fact implies that the time it 
takes to come full circle along the longest loop of size 
C {nRc)~^ is smaller than fs(0 [Eq. (M)] for this 
which yields the condition S/a> (nRcd)^''^ ^ 1. Note 
that this is the same condition at which the effective scat- 



tering time for c? ^ ^ -Rc in Sec. IV A is not renormal- 
ized by the parameter S/a and is given by [cf. Eq. (^0[)]. 
At smaller 6/a, in the interval 1 <C 6/a <^ (nRcd)^^^^ , 
only ^ in the range 1 <^/d< contribute to p^x, 

which gives 



Pxx{B)/pa ^ nRcd ln((5/a) 



(67) 



This equation is valid with decreasing S/a down to (S/a ~ 
1. At still smaller S/a, the particle is scattered out by 
the same AD after passing a distance of order d. Accord- 
ingly, the characteristic time between changes of ADs in 
this regime increases due to the slowing down of the drift 
as Tsa/S. It follows that for S/a 



Pxx{B)/po nRcdS/a 



(68) 



Equations 



3^) match each other at S ^ a'^/d and 
we conclude that the regimes (67), (3^), intermediate be- 
tween those described by Eqs. (22) and (|6^), occur in 
the interval a/d ^ S/a <C (nRcd)^^/'^ . By requiring 
that S at Rc/ls ~ 1 is much smaller than /d we arrive 
at the condition (|64[). The overall behavior of Pxx{B) 
for the case of very weak long-range disorder is obtained 
by adding the "short-scale" contribution, described by 
Eqs. (p2|),(|6^)-(68), and the percolative contribution an- 



alyzed in Sec. |IV D| . This leads to nonmonotonic behavior 
of the MR, such that Pxx{B) first falls off with increasing 
B and then crosses over into the percolative growth. 



V. NUMERICAL SIMULATION 

We have solved numerically the classical equation of 
motion for a charged particle in a random array of hard 
discs in the presence of smooth disorder. In the numerical 
simulation, the latter is characterized by the correlator 



(y(r)y(O)) cx 



1 



[1 + (r/d)2]3/2 



(69) 
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with the ratio d/a 21.5. The concentration n of the 
isolated scatterers has been chosen such that they make 
an important contribution to the resistance at zero mag- 
netic field. Specifically, imcP' ~ 6.2 and l/ls — 0.58, 
where I is the total mean free path and Is — S/Sna. 




FIG. 12. Magnetoresistivity pa;i,(_B) in units of (/i/e^)/fcFZs 
(spin included) as a function of d/Rc for a model system spec- 
ified in the text. Characteristic values of d/Rc are: Rc/ls ~ 1 
at d/Rc ~ 0.25; 5/d = 1 at d/Rc ~ 1.8; nRcd = 1 at 
d/Rc ~ 1.9; nnRl = 1 at d/Rc ~ 2.5; S/2a = 1 at d/Rc ~ 7.1. 

In Fig. |l^ we present the MR data for our model sys- 
tem. Since the characteristic values of the magnetic field 
(given in the figure caption in terms of the parameter 
d/Rc) are rather close to each other, it has not been pos- 
sible to unambiguously separate different regimes. Yet, 
the nonmonotonic behavior of pxx{B), predicted by the 
theoretical analysis, is clearly seen. Note that the size of 
the error bars at d/Rc ^ 5 denotes only the statistical un- 
certainty. In addition to the latter, there is a systematic 
uncertainty originating from a very slow guiding-center 
motion in the large-i3 limit, which makes it difficult to 
observe the true diffusion constant. However, the size of 
the systematic uncertainty is sufficiently small as com- 
pared to the structure of the nonmonotonic dependence 

of Pxx{B). 



VI. SUMMARY 

In summary, we have discussed a rich set of magne- 
totransport phenomena which take place in a random 
ensemble of antidots in the presence of long-range fluctu- 
ations of a random potential. We believe that the model 
studied in the paper adequately describes an antidot ar- 
ray in semiconductor heterostructures with a wide spacer 



for not too high B, when quantum effects (Shubnikov- 
de Haas oscillations) are still weak. We show that even 
weak long-range disorder yields a wealth of pronounced 
effects in the behavior of the magnetoresistance Pxx{B) 
in the antidot array. Essentially, these effects are asso- 
ciated with a magnetic-field induced localization of elec- 
trons which develops when only one type of disorder is 
present. As a result of the localization, Pxx{B) vanishes 
in the limit of large B both in an idealized antidot sys- 
tem without long-range disorder and in a system with 
smooth inhomogeneities without antidots. Conceptually, 
the most striking result of our work is that the interplay 
of two types of disorder does not simply modify the lo- 
calization; in fact, it destroys the localization, so that 
Pxx{B) even diverges in the limit i? — > oo. The diver- 
gence takes place despite a strong falloff of Pxx{B) that 
occurs in intermediate magnetic fields in the case when 
one of the types of disorder is sufficiently strong as com- 
pared to the other. 

Piecing together all the numerous regimes we arrive at 
a rather complex overall picture, due to the interplay of 
several distinctly different mechanisms of the MR. Let us 
list these mechanisms. Some of them are closely related 
to the mechanisms of the MR characteristic to the lim- 
iting cases of strongly non-Gaussian or purely Gaussian 
disorder. Specifically, we have: 

(a) Memory effects operative in the case of strongly 
non-Gaussian disorder (Lorentz gas, or any other sys- 
tem of rare strong scatterers, without long-range inho- 
mogeneities) . These effects lead to a strong negative MR 



(see Sec. HA). In the limit of large B the system is 
insulating. 

(b) Memory effects in smooth Gaussian disorder. In 
this case, the memory effects give rise to a strong posi- 
tive MR, for which, however, the ratio Pxx{B) / po cannot 
be parametrically much larger than 1 , since after having 
reached a maximum value of order 1 it starts to fall off 
with increasing B, thus yielding a strong negative MR 
(see Sec. p^I B| ). In the limit of large B the system is in- 
sulating. 

In the presence of a long-range random potential, the 
mechanism (a) of the MR is destroyed by the diffusive 
motion of electrons scattered by the long-range disorder. 
Nonetheless, as we have demonstrated in the paper, it is 



reincarnated in the form of "hopping rings" (Sec. HIC) 
once the diffusive motion in the long-range disorder turns 
into the drift with increasing B. In this new form, this 
kind of a negative contribution to the MR is developed at 
much larger B as compared to the pure Lorentz-gas sys- 
tem. On the other hand, the mechanism (b) is destroyed 
by scattering on hard scatterers, which checks the nega- 
tive MR associated with t he ad iabaticity of drift in the 
long-range potential (Sec. HI E| ). 

A nontrivial point to notice in our results is that, in 
addition to the above, there are memory effects that are 
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specific to the inhomogeneous system with two types of 
disorder. These are: 

(c) "Diffusion-controlled percolation". As we have 
shown in the paper, scattering by short-range inhomo- 
geneities not just destroys the adiabaticity of motion in 
a smooth random potential, thus checking the strong neg- 
ative MR. In fact, it reverses the sign of the MR by giving 
rise to a positive MR which keeps growing with increasing 
B (Sec. [II B ). This positive MR is a peculiar feature of 
percolation of drift trajectories with superimposed diffu- 
sive dynamics across the drift lines. 

(d) Rcnormalization, by long-range disorder, of the col- 
lision time for hard scatterers. We have demonstrated 
that in a system where the hard scatterers give the main 
contribution to the scattering rate at zero i?, a weak 
smooth disorder can drasti cally suppress the scattering 
rate with increasing B (Sec. IV A ). The effect takes place 
for any type of dynamics of scattering by the long-range 
potential, both for diffusion and drift. The increase of 
the collision time translates into the negative MR. 

The overall behavior of the MR is illustrated in Fig. [T^ . 
Different curves correspond to different n for given Is, II, 
and d. The characteristic field Bad marks the diffusion- 
drift crossover. Figure |l3| describes the case of not too 
weak long-range disorder; specifically, it is assumed that 
Rc/ls <C 1 at i? Bad, which means lL/d<^ (Is/d)^. 

Having listed the main mechanisms of the MR, let us 
now summarize our main results. 

We analyze a "hydrodynamic model" of the chaotic 
antidot array, i.e., tiny antidots scattered with a high 
densityn oo in a smoothly varying random potential 
(Sec. IIIB ). At large B, Pxx{B) turns out to be a grow- 
ing power-law function of the magnetic field. We iden- 
tified several different regimes of the behavior of Pxx{B) 
[Eqs. (|4|)-(|l^)], depending on the parameter Is/V^Il, 
where Is and II are the mean free paths for scattering 
on antidots and long-range disorder with a correlation 
length d, respectively. In the limit _B — > cx3, the hy- 
drodynamic model universally predicts Pxx{B) (x B^^/^^ 
[Eq. (p7|)]. The physics of this divergence is a percola- 
tion of drifting cyclotron orbits limited by scattering on 
antidots. 

Relaxing the conditions of the hydrodynamic model, 
we calculate Pxx{B) in an antidot array of high but fi- 
nite density n (Sec. Ill C ). We show that diffusing cy- 
clotron orbits exhibit intricate dynamics by sticking for 
a long time to a single antidot. This leads to a power-law 
falloff of Pxx{B) in an intermediate range of B, namely 
Pxx{B) oc B~^ In J5 [Eq. (|||)]. The small parameter that 
governs the physics of this transport regime is nRcd, 
where Rc is the cyclotron radius. With further increasing 
B, this mechanism of diffusion is switched off abruptly, 
in a critical manner [Eq. (p3|)]. 




FIG. 13. Schematic behavior of the magnetoresistivity 
pxx{B) on a log-log scale for different values of the concen- 



tration of antidots n: n*^*' > rS"^ > 



> n 



(v) 



keeping 



all other parameters {ls,lL,d) fixed. Only one characteris- 
tic field Bad is shown, at which the crossover between dif- 
fusive dynamics and adiabatic drift in the long-range poten- 
tial takes place. Different curves illustrate different mecha- 
nisms of the magnetoresistance: (i) the magnetoresistance is 
positive owing to the "diffusion-controlled percolation"; (ii) 
due to the adiabatic localization, the concentration of con- 
ducting electrons decreases as In B before the percolation 
becomes effective, which yields a negative magnetoresistance 
Pxx{B) oc _B~^ln_B for intermediate B; (iii) an exponentially 
sharp falloff of pxx{B) at i? ~ Bad (shown as a vertical jump) 
separates the diffusive and drift regimes; (iv) because of the 
memory effects, the collision time for scattering by antidots 
is increased as compared to the Drude value already in the 
diffusive regime {B <^ Bad), which leads to the negative mag- 
netoresistance Pxx{B) oc _B~* for small B; (v) for intermedi- 
ate B, the scattering on antidots stops playing any role and 
Pxx{B) is saturated at a value determined by the long-range 
disorder only, whereas at larger fields the diffusion-controlled 
percolation gives rise to a positive magnetoresistance. 

If the long-range disorder is not too weak, we show that 
the Drude regime and the J5~^ Ini? falloff are connected 
via an exponentially fast decrease of Pxx{B) in a narrow 
range of the magnetic field, which is a trace of the adia- 
batic localization in smooth disorder (Sec. HIE). In this 
regime, Pxx{B) is determined by the interplay of nonadi- 
abatic transitions and scattering on antidots [Eq. (psf)]. 

We calculate the scattering time Tg{B) between col- 
lisions with antidots for extended electron trajectories 
(Sec. |V^, see Eqs. (|^ , (|^) , (|3|) , (|3|) . At large B, r'g 
becomes longer than the Drude time. The smaller n, the 
earlier starts to grow with increasing _B. In a very di- 
lute AD array the rcnormalization of r'g starts already in 
the diffusive regime and leads to a falloff of Pxx(B) cx B^* 
[Eq. (H)]. 
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We solve a scattering problem for a drifting cyclotron 
orbit colliding with an antidot at small 5/a (Sec. iVB). 
We study dynamics of cyclotron orbits skipping along the 
surface of the hard disc and show that the skipping goes 
on almost adiabatically, which results in a strong sup- 
pression of transitions between different drift trajectories 
[Eq. (H)]. 

We analyze complex dynamics of skipping cyclotron 
orbits interacting with an antidot when the ratio Rc/d is 
large, which includes multiple breakaways from the anti- 
dot and multiple returns to it, acco mpan ied by drift along 
closed trajectories in between (Sec. IVC ). We discuss the 
accumulated effect of these multiple collisions with a sin- 
gle antidot in terms of the scattering shift Rh{B) of the 
guiding center of the cyclotron orbit after it finally es- 
capes [Eq. (p7|)]. We also calculate the shift in the limit 
of a smoothly varying environment [Eq. (|5^)] . 

We calculate the percolative magnetoresista nce in the 
antidot array in the limit of small 6/a (Sec. IVD ). A 
variety of different regimes supersede each other with in- 
creasing B, all of which are characterized by a power-law 
behavior of Pxx{B) [Eqs. (|60|)-(|63|)]. The asymptotic be- 
havior of Pxx{B) in the limit i? ^ oo is Pxx{B) oc B^/^^ 
[Eq. (^3|)] , which from a practical point of view to all in- 
tents and purposes is indistinguishable from a saturation 
of the magnetoresistance. This behavior is in sharp con- 
trast to the localization that would develop in the antidot 
array in the absence of long-range disorder. 

We discuss dynamics of cyclotron orbits which stick 
to a single antidot for a long time before hopping to an- 
other one (similar to the B~^ Ini? regime in Sec. Ill C) in 
the case of a very weak long-range disorder (Sec. IV E ). 
Taking this limit eventually restores the B~^ behavior of 
Pxx{B) characteristic to the Lorentz gas. 

We present results of numerical simulations (Sec. 0). 
The numerical data qualitatively confirm the predictions 
of the theory. 
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